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Abstract We derive the life time and loss rate for a 
trapped atom that is coupled to fluctuating fields in the 
vicinity of a room-temperature metallic and/or dielec- 
tric surface. Our results indicate a clear predominance 
of near field effects over ordinary blackbody radiation. 
We develop a theoretical framework for both charged 
ions and neutral atoms with and without spin. Loss pro- 
cesses that are due to a transition to an untrapped in- 
ternal state are included. 

PACS: 03.75.-b Matter waves - 32.80.Lg Mechanical ef- 
fects of light on atoms and ions - 03.67.-a Quantum in- 
formation " 05. 40. -a Fluctuation phenomena and noise 



1 Introduction 

Particle traps enjoy great popularity for the prepara- 
tion and manipulation of coherent matter waves. Promi- 
nent applications are the preparation of non-classical 
states of motion of a single ion IT], the realization of 
quantum gates in quasi-one dimensional ion traps H , 
the transfer of atoms through atomic wave guides f^,^ 
^,y, and the preparation of quantum-degenerate gases 
in electromagnetic-solid state hybrid surface traps Q,^- 
In all these applications, in order to truly benefit from 
the quantum mechanical effects, coherence of the mat- 
ter waves and/or their internal degrees of freedom must 
be maintained as long as possible. Yet, with the phys- 
ical components which provide the trapping potential 
being held at room temperatures, the maintenance of 
coherence seems highly non trivial as the temperature 
gradient between components and trap center may well 
exceed lO^K/m. A careful study of the particles' cou- 
pling to the trap physical components, and the ensuing 
heating of the particles is therefore highly desirable. 

In the past, the heating of single particles in small 
traps has been studied by a number of authors ||9|,|l0|, 
|ll|,|l2|,|l3[|Ij]. As these studies were mostly performed in 
the wake of the recent achievements in ion trapping and 
cooling, the focus in these investigations was on charged 



particles and their coupling to the surrounding metal- 
lic surfaces. In fact, before the advent of laser cooling, 
this coupling provided the dominant cooling mechanism 
for an ion cloud, say, as the low-frequency radiation of 
the ions couples quite efficiently to the lossy currents 
in the metallic trap components M. Yet with the ad- 
vent of laser cooling, temperatures of a few micro-Kelvin 
can be reached which are clearly below the components' 
temperatures, i.e. the particle-component coupling now 
leads to heating, and the trap ground state acquires a 
finite life time. Similar considerations may also be put 
forward for ultracold neutral atoms trapped in miniatur- 
ized traps though the couplings are different: for para- 
magnetic atoms, e.g., they involve fluctuating magnetic 
rather than electric fields close to the trap components. 

In this paper we derive the life time and loss rate for 
a trapped particle that is coupled to fluctuating flelds 
in the vicinity of a room-temperature metallic and/or 
dielectric surface. The theory will be developed for both 
charged and neutral particles with and without spin, 
and loss processes that are due to a transition to an 
untrapped internal state will be included. A detailed 
derivation of previously published results |13] will also 
be given. 

An essential ingredient of the theory are cross- 
correlation functions for thermal electric and magnetic 
fields in a finite geometry. These functions may be sim- 
plified for our purposes because the relevant field fluctua- 
tion frequencies are much lower than the inverse time for 
light propagation from the trapped particle to the sur- 
face and back. It is hence justified to calculate the fields 
in the quasi-static limit, neglecting retardation effects. 
Differently stated, the particle is subject to near field ra- 
diation leaking out of the macroscopic trap components. 
An important consequence is that the near field fluctu- 
ations are much stronger than those of the well-known 
blackbody radiation. This implies larger than expected 
heating rates, as recently pointed out by Pendry m&. 

The paper is organized as follows: in Sec. ||, the model 
is presented in terms of a master equation. We identify 
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the relevant heating and loss rates. Sec. is devoted to 
trapped ion heating. We give the electric field fluctu- 
ations above a flat metallic surface. In Sec. heating 
and loss of a neutral particle with a magnetic moment is 
studied. The final Sec. ra gives a summary and outlook. 
The appendixes contain technical material that is used 
in the main text. 



2 The model: master equation and transition 
rates 

We present here our model for the particle trap and 
its environment (see fig.0, left part). The model is suffi- 
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Fig. 1 Left: trap in front of a flat surface. Right: heating 
and loss processes. 

ciently simple to allow for analytical calculations of the 
relevant heating and loss rates, but also reflects a typi- 
cal experimental geometry. We consider a single particle 
bound in a harmonic trap potential whose center is lo- 
cated at a distance z from an infinite flat surface. We 
consider that this distance is much larger than the size 
of the particle's center-of-mass wave function. In this 
regime, the overlap with the surface is negligible, and the 
coupling to the surface is mediated via electromagnetic 
fields. We also focus for simplicity on a single degree of 
freedom in the harmonic well. 

The heating of the particle is described by the tran- 
sition rate /q^i from the trap ground state |0) to the 
first excited state |1) (see fig.0, central part). In sub- 



section 2.1, such a 'heating rate' is determined from a 
master equation for the particle's motion in terms of 
harmonic-oscillator matrix elements, on the one hand, 
and the spectral density of a fiuctuating force field, on 
the other. 

As a second application, we investigate loss processes 
in magnetic or optical traps where only a subset of inter- 
nal states is trapped (see fig.|^, right part). This model 
describes magnetic traps, for example, where only low- 
field-seeking Zeeman sublevels can be trapped. A loss 
process occurs when a fluctuating field induces a fiip 
\i) -^ I/) of the particle's internal state. We assume 
that the particle is then rapidly expelled and lost from 
the trap. The relevant loss rate il-t/ is given in sub- 



section 2.2 in terms of internal matrix elements for the 



particle's magnetic moment, on the one hand, and the 
magnetic field fluctuation spectrum, on the other. 



2.1 Heating 

As mentioned before, we focus on the heating of a single 
degree of freedom for the trap vibration. The displace- 
ment X of the particle relative to the trap center r is 
chosen along the unit vector n and written in terms of 
a creation operator b. The interaction potential reads 

V{r, t) = -x- F(r, t) = -a{b + b^) n • F(r, t), (1) 

where a = (?i/(2Af 12))^/^ is the size of the trap ground 
state {M is the particle mass and il the trap frequency) 
and F(r, t) the force acting on the particle. This force is 
fluctuating, and it is convenient to use a reduced density 
matrix description for the particle when the force fluc- 
tuations are averaged over. The density matrix p evolves 
according to a master equation that is written in eq. (E2) 
of appendix |A.1| for a general coupling. For the Hamilto- 
nian (Rl), we get the following relaxation dynamics p7| 
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— {bb^p + pbb^ - 2b'fpb) . 



(2) 



In this equation, the transition rates j± = 7(r; ±f2) are 
proportional to the spectral density Sp of the force fluc- 
tuations taken at the trap vibration frequency fi 



7(r;w) 



^n,njS')?(r;w), 



(3) 



where Sp (r; cu) is defined by 

+00 

S'^{r;io)^ [dT{F,{r,t + r)F,ir,t))e^-\ (4) 



From the master equation (0), it is easy to obtain 
rate equations for the populations of the trap levels. For 
the ground state population poo = (0|/5|0), we get 



POolrclax = -7-POO + 7+Pll- 



(5) 



Note that the transitions towards higher (lower) trap 
levels occur with a rate equal to 7_ (to 7+). In particular, 
the quantity 7_ gives the depletion rate of the ground 
state population. The heating rate we are interested in 
thus equals 



/l)^i(r) =7- = -2 Vn,nj5]^(r;-/2). 

h 



(6) 



Note that the same result may be obtained from Fermi's 
Golden Rule, by assuming a mixture of initial states for 
the fluctuating force field and summing over its final 
states. In Sees. || and ^, the heating rates for trapped 
ions and spins are computed using (jm. The main goal 
of the calculation is therefore the spectral density of the 
relevant force (electric or magnetic fields). 
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Finally, the master equation (0) also allows to de- 
scribe the decay of the coherences between trap states 
which is a hazardeous process for quantum bit manip- 
ulations. The coherence between the lowest trap levels 
relaxes according to 
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7+ +7- 



relax 



POl + V27+P12 



(7) 



the metallic surface nearby. The force in the interaction 
Hamiltonian (nh is given by the electric field 



F(r,t) = r7E(r,i) 



(10) 



where q is the ion's charge and r the position of the trap 
center. 



We see that the coherences decay with a similar rate as 
the populations. This is a consequence of the interaction 
Hamiltonian (]^, and different results are obtained using 
other couplings or adding explicit phase noise, see, e.g., 
Refs. |l^ , [l3|| . In the following, we focus on the population 
dynamics for simplicity. 



2.2 Internal state 



In magnetic or optical traps for neutral particles, the 
trap potential depends on the internal atomic state (see 
fig||, right part). If this state is changed due to fluctua- 
tions in the magnetic field, the particle may be subject to 
an anti-trapping potential and strongly perturbed. The 
interaction Hamiltonian for spin flips |z) — > |/) is the 
Zeeman interaction 



Vz{r,t) = -fi-B{r,t), 



(8) 



where fi is the particle's magnetic moment and B(r, i) 
the fluctuating part of the magnetic field. For this inter- 
action, a master equation similar to (H) may be formu- 



lated from the general theory outlined in appendix A.l 



This equation is not very instructive, however, if we as- 
sume that the particle is lost as soon as it reaches the 
state I/). In this case, it is sufficient to quote the tran- 
sition rate il_>/ obtained from ( |42| ) 

-f^-'/W=V 72 Sg'^iri-ujf,), (9) 



where Sg is the magnetic field fluctuation spectrum de- 
fined by an expression similar to (0) , and fiujfi = Ef — Ei 
the energy difference between initial and final internal 
states. (We switch to greek subscripts to avoid confu- 
sion with the initial state label.) In a magnetic trap, e.g., 
\i), I/) are magnetic sublevels and the frequency tufi a 
Larmor frequency in the bias field of the trap. In opti- 
cal traps, we consider the hyperfine components of the 
atomic ground state, LOfi is thus the hyperfine splitting. 



3. 1 Electric field fluctuations 

In the formula (ph for the heating rate, we need the spec- 
tral density of the electric field fluctuations ^^(r;^). 
This quantity is conveniently obtained by making use 
of the fluctuation-dissipation theorem outlined in ap- 



pendix A. 2. According to this theorem, the field's spec- 



tral density is proportional to the imaginary part of the 
field's Green function Gij(r,r;a;), multiplied with the 
Bose- Einstein mean occupation number (eq.(|4q)). The 
geometry we have chosen is sufficiently simply to allow 
the Green function to be calculated analytically ||l^ . Re- 
call that the Green function describes the electric field 
radiated by an oscillating dipole (cf. eq.(|47|)). This field 
is the sum of the dipole field in free space plus the field 
reflected from the surface. The free space field leads to a 
term G^ (r, r; w) in the Green function that is actually 
independent of the trap position r; it gives the spectral 
density of the blackbody field (the Planck law) 
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(11) 
(12) 



where T is the temperature of the surface (we put the 
Boltzmann constant ks = 1). 

To calculate the field reflected from the surface, we 
expand the free space dipole field in plane waves and 
apply the Fresnel refiection coefficients rs.p(u) for each 
wave incident on the surface (s and p label the two trans- 
verse field polarizations and u is the sine of the angle of 
incidence). The resulting Green function GJ" (r,r;a;) 
characterizes the modification of the thermal radiation 
in the near field of the surface. The radiation density 
is increased with respect to the far field expression ( pl| ) 
because it also contains non-propagating (evanescent) 
waves. The corresponding spectral density depends only 
on the distance z to the surface and may be written in 
the form fll 



3 Heating of a trapped charge 

In this section, the master equation of the previous sec- 
tion is applied to the most simple situation, that of an 
electrically charged particle in a harmonic trap |9|,[iO|,PT|, 
|lj,|l3|,|l4|. As mentioned in the introduction, the ion is 
heated up because fluctuating electric fields leak out of 



S 



{nf)ij 



{r;uj) = S^E^\uj)gij(kz) 



(13) 



where the diagonal tensor g^ has the dimensionless ele- 
ments gxx = gyy = g\\ and g^^ = 5^ with {k = \uj\/c) 



4-oc 



5||(A:z) = I Re / ^e^^'^^ {r,{u) + [u^ - l)r,(^)) 
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+00 
3 f u^du 

9±[kz) = 2 / ~^' 



2ikz 



'rp{u), 



vT^^??, < u < 1, 

iy/u^ — 1, M > 1. 



Finally, the relevant Fresnel coefficients are 



(14) 
(15) 



rp{u) = 
rs{u) = 
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V + \/e — u^ 



(16) 



where £(aj) is the relative dielectric function of the bulk 
metal. 

For typical trap frequencies the corresponding elec- 
tromagnetic wavelength is much larger than z, so we 
can restrict our calculations to the quasi-static limit 
z <^ X and find analytical expressions for the tensor 
elements (|l^). The details are outlined in appendix p. 
We have to distinguish between the case of a large and 
a small skin depth of the conducting material compared 
to the distance z. The skin depth, which is the charac- 
teristic length scale on which an electromagnetic wave 
entering a conducting solid is damped, is given by (for 
OO) |I| 

5 =^^2^^ (17) 

where g is the specific resistance. Since in our frequency 
regime the dielectric function for a metal is dominated 
by the zero-frequency pole, it is related to the skin depth 

by 

eH « ^ = tSi- (18) 

In appendix B.l, we derive approximations for the 
functions g\\,±^{kz) in the form of inverse power laws 
(eqs.( |5^ , |54| )). Both regimes of large and small skin depth 
can be covered by the following interpolation formula 



5y {kz) 



3J2 



'ij + Oij 



(19) 



where Sy is a diagonal tensor with the elements Sxx = 
Syy = i, Szz — 1- Thus we arrive at a final expression for 
the electric field spectrum, applying the high tempera- 
ture limit of the Planck law (|lj): 
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(20) 



We note that in the case of a short distance, the parallel 
and perpendicular tensor elements both show a 1/z^- 
dependence and differ by a factor of 2, whereas for larger 
distances the tensor elements are equal and show a l/z^- 
behavior. 

The 1/z'^ power law of the regime z <C (5 may be 
understood in terms of image theory: the electrostatic 
dipole field varies precisely as 1/r^ and its refiection from 



the surface is characterized by the factor (e— l)/(e+ 1) w 
1 + i{k5)'^. The imaginary part of the reflected field thus 
reproduces (p^). This is the regime discussed in Ref . |p5[ . 
It is interesting to note that for a larger distance z 3> 
S, the field fluctuations are enhanced with respect to 
the electrostatic regime (see fig. pi). This is due to the 
fact that the dipole field is more efficiently damped in 
the conductor because the exponential decay in the skin 
layer quenches the algebraic penetration of the field. 

For completeness, we also mention the limiting case 
of a perfectly conducting surface (e — > 00) whose skin 
depth S vanishes. The previous asymptotic expansion 
does not cover this case. The coefficients g\^±{kz) given 
in the appendix H, eq.(pq), show damped oscillations 
with a period equal to the wavelength. In the short- 
distance limit z ^ A, we get g\\{kz) -^ — 1 and g±{kz) -^ 
1, the divergence at z — > thus disappears. The electric 
field fluctuations are essentially those of the free space 
blackbody spectrum, with a minor modification due to 
the boundary conditions. 



3.2 Heating rate 

We plot in fig. the heating rate <M) for an ion (trap 
frequency J7/27r — 1 MHz) above a copper surface. The 
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Fig. 2 Heating rate for a trapped ion. Dots: coupling to 
electric proximity fields, computed from (|l5|). The solid line 
is obtained using the asymptotic formula (H). Dashed line: 
coupling to thermal voltage fluctuations. 
Parameters: trap frequency J7/27r — 1 MHz, copper substrate 
with g = 1.7 X 10^^^ cm at T = 300 K. The ion mass is 
M = 40amu, and its charge q = e. The trap axis is perpen- 
dicular to the surface, n = e^. The thermal voltage fluctua- 
tions are characterized by a circuit resistance 1 Q, lid] . The 
endcaps are separated by twice the ion-surface distance. Size 
and inverse lifetimes of typical ion traps are indicated by the 
shaded rectangle [0,^|2J. 



dots are based on an exact (numerical) evaluation of the 
^-coefficients ([l5|), while the solid line uses the interpola- 
tion (19). The change in the power law at the skin depth 
is clearly visible. Note the marked increase of the field 
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fluctuations compared to the free space blackbody level 
(dotted line). Also shown is the estimate given by Lam- 
oreaux pQ l who modeled the trap in terms of a resistively 
damped capacitor with a thermally fluctuating voltage 
(Johnson noise). Wineland et al. Cj] pointed out that 
realistic estimates for the corresponding resistance actu- 
ally give smaller heating rates. Our results suggest that 
the miniaturization of ion traps down to fiia sizes en- 
tails difflculties to maintain long coherent storage times, 
unless all physical components are cooled down. 



4 Trapped spin coupling to magnetic fields 

In this section, we turn to traps for neutral particles and 
consider the Zeeman coupling m) of the atomic magnetic 
moment to a fluctuating magnetic field. In magnetic and 
optical traps, this coupling may induce a spin flip to a 
non-trapped state (magnetic sublevel or hyperfine state) . 
This implies a nonzero los s rate from the trap that we 
calculate in subsection 4.1. On the other hand, the Zee- 



man interaction also exerts a force proportional to the 
gradient of the magnetic field. If this force fiuctuates, it 
does not necessarily flip the atomic spin, but excites the 
atom into a higher trap level. The corresponding heating 
rate is the subject of subsection 4.2. 



4-.1 Spin flips 

4.1.1 Magnetic field correlations. We first compute the 
magnetic field fluctuations in the vicinity of the solid sur- 
face. By analogy to the ion case, we use the fluctuation- 
dissipation theorem ( [48|) and determine the Green tensor 
for the magnetic field. In fact, the calculation is very sim- 
ilar to that for the electric field: starting from the field 
radiated in free space, we expand it in spatial Fourier 
components and compute for each plane wave the reflec- 
tion at the solid surface. It turns out that the Fresnel 
coefficients for the magnetic field are identical to those 
for electric fields, except that one has to exchange the s- 
and p-polarizations. We thus get the following near-field 
correction to the magnetic field fluctuation spectrum 
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(21) 



Similar to (|l^), hij is a dimensionless and diagonal ten- 
sor with elements 

mkz) = ^Re I !i^e2''=- {rp{u) + {u' - l)r,(u)) , 



+ OC 



h^{kz) = -Re / !i^e=^"=-V,(u). 
2 J V 



(22) 



For experimentally relevant parameters, the magnetic 
fields at the resonance frequency have a wavelength 



(at least some cm) much longer than the size of the 
trap. This implies again that we need the short-distance 
asymptotics 2; ^ A of (pTj). A calculation outlined in 



appendix B.2 gives the following interpolation formula 



that covers both regimes of a large and small skin depth 
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(23) 



where Sy is the diagonal tensor introduced in (pj|). The 
magnetic field spectrum ( |2l] ) thus equals in the high- 
temperature limit 
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(24) 



Note the different exponents for the distance dependence 
compared to the electric field fiuctuations (po|). 

If the trap distance is small compared to the skin 
depth, we recover the magnetic field spectrum given in 
eq.(lO) of ||l^, apart from the fact that the parallel ten- 
sor components (^s^x-, Syy) differ. This difference is due to 
the fact that the calculation of [|5| uses the Biot-Savart 
law to get the magnetic field from a statistical model of 
polarization currents in the solid. This approach is valid 
for stationary currents only, and a difficulty appears at 
the surface because the model for the currents is not 
divergence-free there. Therefore, while the magnetic field 
perpendicular to the surface is correctly described, the 
parallel components are overestimated. 

4-. 1.2 Internal matrix elements. In order to compute 
the spin flip loss rate we have to evaluate matrix ele- 
ments of the total magnetic moment operator as indi- 
cated in (0). This operator is in general given by 



A*^ 



/is .9lL + gsS - gi 1 , 



(25) 



with fig the Bohr magneton, L the total orbital angular 
momentum operator, S the electronic spin operator, I 
the nuclear spin operator and gL, gs and gi the corre- 
sponding (/-factors. Since the proton mass rup is larger 
than the electron mass m^ by three orders of magnitude, 
we can neglect the contribution of the nuclear magnetic 
moment. Furthermore, the reasonable restriction to an 
atomic ground state with L = reduces the problem to 
the calculation of matrix elements of solely the spin op- 
erator. Together with the fact that the tensor hij in ( |2^ ) 
for the magnetic field correlations is diagonal, we can 
focus on terms of the form 



K/lA^a 



^^Ws\{f\So.\^)\' 



(26) 



In the following we will restrict ourselves to two ex- 
treme cases: the coupling between two Zeeman sublevels 
in the presence of an external magnetic field and the 
coupling between two hyperfine ground states without 
external fields appfied. The former case is e.g. reafized 
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in a magnetic trap, whereas the latter corresponds to 
optical traps. 

In the case of a magnetic trap the trapped atom is 
subject to a constant magnetic field with strength Bq in 
the center of the trap, assuming the atom is not moving. 
The magnetic sublevels are split due to the Zeeman effect 
by the Larmor frequency lol = gsl^BBo/h. (We focus on 
a vanishing nuclear spin for simplicity.) Without loss of 
generality we can assume the magnetic field to be lying 
within the xz-plane, since the diagonal tensor in ( p3| ) has 
the symmetry property hxx = hyy . If the magnetic field 
forms an angle 9 with respect to the z-axis, we denote 
by \m)e the basis states with quantization axis parallel 
to the magnetic field (the 'trap basis'). Rewriting ( pq ) 
leaves us to calculate matrix elements of the form 



l(/lMaN)l^ ^ fJ'l9l\e{mf\Sa\m^)g\ 



(27) 



These elements are evaluated by expanding the spin vec- 
tor components in a rotated coordinate system (denoted 
by the prime) adapted to the trap basis. The result is 
the following: 



e{mf\Sx\m.i)g = (g{mf\S'_^\mi)e + e{mf\S'_\m.i)g) 

+ e{mf\S'^\mi)e sin 9, 
e{'mf\Sy\m.,)g = - (g{mf\SL\mi)e - g{mf\S'^\mi) 

e{mf\S3\m.i)g = (g{mf\S'^\mi)0 + g{mf\S'_\mi)g) 



cos 9 



j(to/|S'3 
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3l"H 



2 

(28) 



where Sg is the z-component of the spin operator and 
S'_^, S'_ correspond to raising resp. lowering operators in 
the trap basis, whose action is known |2^. In the case 
of an electronic spin S = 1/2, the trapped (untrapped) 
level is the \Tni)g ~ \~l/2)g {\mf)g = |l/2)0) Zeeman 
sublevel, respectively. The matrix elements (|8|) then be- 
come 



(l/2|5,|-l/2)e 
(l/2|5,|-l/2>, 
(l/2|5,|-l/2)e 



cos 9 

^~' 
i 

^2' 
sin6' 



(29) 



With this result, we can compute the magnetic loss 
rate (34) below. 

In the case of an optical trap we have to take into 
account that the nuclear spin couples to the electronic 
spin, F = S -I- I, and causes the ground state to spht 
into hyperfine levels, separated by a frequency ujhf- We 
are now interested in the transition probability from one 
hyperfine ground state to another. Thus, for this case we 
can write (Effl as 



(30) 



\{f\^^c.\^)\'^f^Wsm\S^\F,)\' 



A transition from one hyperfine ground state to an- 
other can take place between different magnetic sub- 
levels. Thus we first have to calculate the transition rate 
between two of these states. This is done by expand- 
ing the basis states in the uncoupled basis, choosing the 
quantization axis taken along the z-axis: 
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(31) 



where C^^' are the Clebsch-Gordan coefficients. The 
matrix element between two hyperfine magnetic levels is 
then 



{Ffmf\S^\E,7 
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^Ffnif ^FiTTii 



{m's\S^\ms). (32) 
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Note that the nuclear spin does not flip in the transi- 
tion. Again the action of Sa onto the electronic spin 
states \ms) is well-known in (|3^). We obtain an effec- 
tive transition rate between the two hyperfine manifolds 
by summing the rates over all final 77i/-levels and tak- 
ing the average over the initial mj-levels. This gives the 
following result for the hyperfine matrix element (30) 



\{Ff\S^\F^ 



1 



2F^ + l 



E \{Ffmf\Sa\F^m, 



(33) 
We finally note that this calculation assumes that the 
frequencies for the transitions \Fimi) — + \Ffmf) are all 
equal to the hyperfine splitting ujhf ■ This is a good ap- 
proximation if Tildhf is large compared to the optical 
trap potential (that may lift the degeneracy of the hy- 
perfine states even without a static magnetic field). 



^. 1.3 Loss rate. Combining the matrix elements (E6 
for the magnetic moment, the magnetic field spec- 
trum (p4) and eq.(Bh, we get the following loss rate for 
a magnetic trap 

^-/ W = ^¥^$Y J2(^Ukz) + 1) \{f\S^\^)f ■ 

(34) 
For the case of an electronic spin 5* = 1/2 and no nuclear 
spin we can use the matrix elements from (G3) and obtain 



r_ 



.« 



i^Ws^lT 



Uneofi^c^ 



{(ft||(fcz) + l)(l- 
+ {h^{kz) + l)sm'^9} 



cos 



(35) 



This loss rate is plotted in fig. y| for two different Larmor 
frequencies lul , with the trap bias field chosen parallel to 
the surface {9 = 7r/2). We see that quite large loss rates 
occur if the trap center approaches the surface down to a 
few micrometers. Again, miniaturized traps have to face 
the influence of larger noise flelds. 
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Fig. 3 Loss rates in a magnetic trap above a copper surface. 
Dots (solid lines) : results based on (021) (on the asymptotic 
interpolation (p3|) ) . Results for two different Larmor frequen- 
cies u}l/'2it = 1 MHz (curve a) and 100 MHz (curve b) are 
shown. The arrows mark the corresponding skin depths. The 
shaded area indicates experimental data obtained in Kon- 
stanz and Heidelberg MM. 

Parameters: spin S — 1/2, magnetic bias field aligned par- 
allel to the surface. The loss rate due to the blackbody field 
(the prefactor in (H)) is about 10~"s~^ at 100 MHz (not 
shown) . 



In fig. H, we plot the loss rates obtained from the ef- 
fective matrix element ( p3[ ) for hyperfine-changing tran- 
sitions. The data are calculated for the lower ground 
states of trapped ^^Rb and -'^^'^Cs. One observes that 
these rates are much smaller than those for magnetic 
traps. It is interesting that this reduction is due to the 
skin effect: indeed, the magnetic field fluctuations (p3) 
in the intermediate-distance regime 5 <C z <C A are pro- 
portional to (5^ ex w"'^/^. Larger transition frequencies 
thus lead to smaller loss rates. 
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Fig. 4 Loss rates due to hyperfine-changing transitions in 
an optical trap above a copper surface. Dots (solid lines): 
results based on (e3) (on the asymptotic interpolation (E3|)). 
Results for two different atoms are shown: ®^Rb (7 = 5/2, 
ujhf/2-k = 3.04 GHz, transition Fi = 2 ^ 3 ^ Ff) and ^^^Cs 
(/ = 7/2, ujHF/2n = 9.193 GHz, transition F, = 3 ~* 4 = 
Ff). The horizontal dotted line marks the corresponding loss 
rates in the free space blackbody field. 



4-.2 Heating of the cm. motion 

This case is treated by analogy to the trapped ion. 
The Zeeman interaction (g) gives the following magnetic 
force 

Fz(r,t)=V(/x.B(r,t)) (36) 

that couples to the displacement of the particle from its 
equilibrium position. The matrix elements for the dis- 
placement are that of a ID harmonic oscillator and are 
given in subsection 2.1. We are left with the calculation 
of the magnetic force's spectral density. To this end, re- 
call the identity 



{Fz^ir,t')Fz,{r,t)) 
d d 



dru dr2j 



{Vzirut')Vz{r2,t')) 



(37) 



The relevant information is thus contained in the cross 
correlation function for the magnetic field at two differ- 
ent positions ri 2. From the fluctuation-dissipation theo- 



rem (appendix A_^), this correlation function is propor- 
tional to the Green function Hij {ri,r2;ui) for the mag- 
netic field. To simplify the calculation, we focus on a 
trap with an axis n perpendicular to the surface. Ac- 
cording to (|g), we then only need the zz-component of 
the force fluctuation tensor. In the identity (|37|), it is 
thus sufficient to take two positions ri^2 = (R-, ^1,2) that 
differ only in the vertical coordinate (R = {x,y) denotes 
the coordinates parallel to the surface). It may now be 
shown that the surface-dependent part H^J' {ri,r2;cj) 
of the Green tensor depends only on the average distance 
z = (zi -|- Z2)/2 and the lateral separation R2 — Ri p8[ . 
This is clear, e.g., from image theory. Since Ri = R2 for 
our special case, we may write 



Hlf\R,zi,R,z2;ij)^H, 



(«/), 



R, z, R, 



z:uj 



(38) 



where the right-hand side is the Green function taken 
at identical positions that has been calculated in subsec- 
tion |4.1.1| . 

We now use the results (p6, 57) for the magnetic cor- 
relation tensor (app. B.2), write z = {zi + Z2)/2 and 
differentiate with respect to zi^2- All told, both asymp- 
totic regimes of small and large skin depth are described 
by the interpolation formula 



SFA^;io) 



fi^T (z|m2 



A^iN) 



647rp 



1 



1553 



(39) 



This spectrum is already summed over all flnal Zeeman 
states, assuming that all of them are trapped. The aver- 
age for the magnetic moment is taken in the initial state. 
For an atom with L = 0, S = 1/2 in the ground state, it 
equals 5|/Xg ~ 4/ig where ^b is the Bohr magneton. 

If the trap distance is small compared to the skin 
depth, we recover the expression (11) of ||l^ for the heat- 
ing rate 



ro^i(r) 



GAnhnMgz^' 



(40) 
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apart from different weights for the paraUel and per- 
pendicular spin components. This is due to the different 
magnetic field correlation tensor (E4h that has already 
been discussed above. 

In fig. |l we plot the heating rate Iq^i obtained from 
the magnetic fluctuation spectrum (p9) for a typical trap 
above both a copper and a glass surface. The heating 
rate above glass is much smaller because glass is a poor 
conductor. For a copper substrate, note the crossover 
when the distance becomes larger than the skin depth. 
A remarkable result is the large value of the heating rate 
for small traps (dimensions below the yum range). 
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A Statistical tools 



10' s" 
1 s" 



tn 10-*S" 



?P in-12.-l 



10"'-s" 



10"''*s" 







' : 


^^^^M 




^^^--^ 


Typical micro - traps \ 


Spin 


heating (Cu) 


^- : 


^ Spin 


heating ( glass ) 


skin depth (Cu) 



100 nm 



1 (im 10 (im 100 |im 

distance from surface 



1 mm 



Fig. 5 Heating rate for a trapped spin above copper and 
glass substrates. 

Parameters: trap frequency i7/27r — 100 kHz, M = 40amu, 
magnetic moment fi = hb = 1 Bohr magneton, spin S = 1/2. 
The heating rate due to the magnetic blackbody field (not 
shown) is about 10~^®s~^. For the glass substrate, a di- 
electric constant with Re £ = 5 and a specific resistance 
g — lO^'^ncm are taken. These values are used in the short- 
distance asymptotics (M) to compute the magnetic field fluc- 
tuations. 



5 Summary and outlook 

To summarize, we have developed a theoretical frame- 
work for the systematic investigation of the heating and 
concomitant loss of coherence in small particle traps. 
Our results indicate a clear predominance of near field 
effects over ordinary (free space) blackbody radiation. 
They establish upper bounds for life times in a variety 
of experimentally relevant types of traps. 

The present model is restricted to particle motion 
in a single dimension, and the extension to a three- 
dimensional trap geometry is an obvious step for future 
work. A theory beyond the rate equations discussed here 
could include noise-induced shifts of the particle's energy 
levels. Finally, still other interactions might be consid- 
ered for neutral atoms. The coupling to electric fields via 
the polarizability tensor is currently under investigation. 



A.l Master equations 

We outline here a general master equation n^ that de- 
scribes the reduced dynamics of a system coupled to a 
reservoir. The coupling Hamiltonian is given in terms 
of an arbitrary system operator s, a fluctuating force 
F(r, t), and a coupling constant g 



V{r,t) = -gs-F{r,t). 



(41) 



Throughout this paper, the parameter r denotes the trap 
center position. For a trapped ion, e.g., the system opera- 
tor s would describe the displacement of the ion from the 
trap center, see eq.(H). In the Markov limit and ignoring 
reservoir-induced level shifts, the relaxation dynamics of 
the reduced system density matrix p is 



P|rc 



g^^ S'^{r;Lo) 



-24+V4- 



S'^ir; 



1-3 



-24-V^^O 



^^t'p- 



'Jf'^ 



p^'^ 



p^^-^ 



(42) 



where the s*^^) is the positive (negative) frequency part 
of the system operator. More precisely, the free system 
evolution in the Heisenberg picture is given by 



5(i) = s(+)e-*'^*-Fs(-)e*'^* 



(43) 



where huj (> 0) is the energy difference between two ad- 
jacent system states. The spectral density in (E2h is de- 
fined by (cf. eq.(|)) 



+ 00 

S%{j-L^)= JdT{F,ir,t + T)F,ir,t)) e'^\ 



(44) 



The master equation (42) allows to derive rate equations 
similar to (||), and these show that the rates proportional 
to Sp (r; +uj) govern spontaneous and stimulated decay 
processes, while excitation processes are proportional to 
Sp{r;—uj). The latter correlation function is thus rele- 
vant for our heating problem. 



Loss and heating of particles in small and noisy traps 



A. 2 Fluctuation-dissipation theorem 

In a reservoir at thermal equilibrium, there is a relation 
between the cross correlation tensor for the field fiuc- 
tuations and the field's Green tensor iQ. This relation 
also holds for correlations taken at different positions in 
space, that we have to compute in subsection 4.2. For a 
force field F(r, i), the cross correlation tensor is defined 
by generalizing ( [44| ) 



+00 
^)?'(ri,r2;w)= / dr (F,(ri, i + T)F,(r2, t)) 



(45) 



The Green function is defined as the force field created 
by a classical monochromatic, localized disturbance a at 
ro (e.g. the electric field of an oscillating point dipole). 
The interaction Hamiltonian density is 



e-"^'6{r-ro)a-F{r,t). 



(46) 



In thermal equilibrium, the average linear response to 
this source is a harmonic field (F(r,t;ro)) that depends 
parametrically on the source position ro and is propor- 
tional to the displacement a. The Green function is the 
corresponding proportionality factor 



(F,(r,t;ro))=e-''^*^G,,(r,ro;c^)a, 



(47) 



(The averaging (• • •) removes the oscillations of the free 
field.) The fluctuation-dissipation theorem now states 



S'-^iri,r2;u) = ^— — ^-^.... ^^ 



ImGy(ri,r2;w). (48) 



Note that in terms of the mean thermal occupation num- 
ber nth = l/(e^'^/^ - 1), one has (for uj > 0) 

S'^iri,r2;Lu) = 2h{nth + 1) Im G,,(ri, r2;a;), (49) 



S'/{ri,r2;-uj) = 2^nthlm Gij{ri,r2;uj). 



(50) 



At zero temperature, fith — 0, and only the first line 
survives. The relaxation dynamics is then entirely due to 
spontaneous decay, induced by the vacuum fluctuations 
of the force field. Heating processes are suppressed. At 
high temperature, fith ^ 1, the fluctuation spectrum 
becomes independent of the sign of w. In the master 
equation, decay and excitation rates are then nearly the 
same. 



B Asymptotic expansion of electromagnetic 
field spectra 

B. 1 Electric field 

We outline here the asymptotic expansion for the coef- 
ficients 911 x(fcz) that characterize the electric field fluc- 
tuations ( |l3| ) in the near field kz <^\ oi the surface. 



The inspection of the integrals (E5[) shows that the 
exponential e^*'^^" « g-2fczti decreases on a large scale 
u ~ l/{kz) ^ 1. On the other hand, the other factors in 
the integrands increase as powers of u. The value of the 
integral is thus dominated by values u ~ Wmax around 
the maximum Wmax ~ ^/(kz) 2> 1. It is therefore accu- 
rate to use asymptotic expansions of the Fresnel coeffi- 
cients for large u 3> 1. The asymptotic form of the co- 
efficients depends, however, on whether u'^^y. is smaller 
or larger than the magnitude |e| of the dielectric con- 
stant. These two regimes are discussed in the following. 
Their physical significance follows from the relation ( |l^ ) 
between s and the skin depth S. 

The limit 1 <^ |e|^^^ ^ u corresponds to a distance 
small compared to the skin depth, z <C (5 ^ A. In this 
regime, we get the following asymptotic expressions for 
the Fresnel coefficients dlq) 



rp{u) 
rs{u) 



e-1 

FTT' 

£-1 

4u2 



(51) 



The integrals (IT^) for the tensor elements are then eval- 
uated to 



^"(^")"l6(fcz)^ 
g±{kz) « 2g\i(kz). 



3 , e- 1 
im 



3,52 



e + 1 16kz^ ' 



(52) 



In the opposite limit of a small skin depth, i.e. (5 <C 2; <C 
A, we have 1 <C u ^ kl^^^i a-nd the reflection coefficients 
show the asymptotic behavior 



rp{u) -^ 1 + 



2i 



ude 



rs{u)^ -1 + ^. 

This yields tensor elements of the form 

3 „ 1 3(5 



(53) 



5||(fcz) ; 
g±{kz) 



— Re 

4(/fcz)2 ^^ ^ 

'9\\{kz). 



8kz 



(54) 



The regimes (|5^ , |54|) a re readily combined into the inter- 
polation formula (|l9|). 

In the limit of a perfectly conducting (pc) surface 
(e — > c»), the skin depth S vanishes, and the refiection 
coefficients^ 16|) are equal to rp^s ~ ±1 (cf. eq.(p3|)). The 
integrals (|l5) may be evaluated explicitly, and one gets 



PC 



9\\ikz) 
g±{kz) 



3 /sin2fcz cos2fcz sin2fci 



2 V (2fcz)3 
' sai2kz 



{2kzY 
cos2kz 
(2fcz)3 ~ (2fcz)2 



2kz 



(55) 



Note that these functions have finite limiting values at 
z ^ 0, which is different from the behavior (52) above a 
surface with a finite conductivity. 
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B.2 Magnetic field 

The asymptotic evaluation of the coefficients hu^{kz) 
for the magnetic field spectrum (K^) proceeds similar to 
the case of the electric field. 

For a skin depth larger than the trap distance, we 
expand the reflection coefficients in the regime 1 <C 
|e|"'^/^ <C u. The asymptotics of the tensor elements ( p2[ ) 
is then given by 



z< (5< A 

h\i{kz) 

h±{kz) 



3 ^^^ (e-f)(£ + 5) 



32kz 
3 



16k^S'^z ' 



16kz 



Im(e- f) « 2h\\{kz) 



(56) 



We used the approximation |e| 3> f appropriate for a 
good conductor. 

In the opposite limit of a small skin depth, we find 



(5< z< A 

hu{kz) 



9 ^ 

I6(fcz)4 '^'^ ^ 



96 



h±{kz) w 2/i||(fcz). 



(57) 



Both expressio ns (|5^ , p7| ) are reproduced by the interpo- 
lation formula (^31). 
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